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Using Krichever-Phong's universal formula, we show that a multiplicative repre- 
sentation linearizes Sklyanin quadratic brackets for a multi-pole Lax function with 
a spectral parameter. The spectral parameter can be either rational or elliptic. As a 
by-product, we obtain an extension of a Sklyanin algebra in the elliptic case. Krichever- 
Phong's formula provides a hierarchy of symplectic structures, and we show that there 
exists a non-trivial cubic bracket in Sklyanin's case. 
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A starting point for many soliton systems is a Lax equation 

L=[P,L], (1) 

where L and P are operators. 

For most finite-dimensional integrable systems, this equation can be interpreted as a flow 
on the space £ of meromorphic matrix functions L(z) on some Riemann surface V (typically, 
a Riemann sphere or an elliptic curve), where the positions of the poles are fixed. 

The algebraic-geometric procedure for constructing the exact solutions (see |9J for a brief 
outline and additional references) is merely a parametrization of the space £ in terms of a 
spectral curve with marked points and a divisor. 

A spectral curve is defined by the equation 

f : R(k, z) = det (L(z) - kl) = 0. (2) 

For L(z) in general position, it is a smooth Riemann surface of finite genus. For every 
point Q = (k, z) of T, there exists the unique eigenvector ^{Q) of L(z) satisfying 

L{z)i>(Q) = kiP(Q), 
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normalized so that its first component is one, i.e., ipi(Q) = 1. The vector function ip(Q) is 
meromorphic on T. Due to the Lax equation, the spectral curve does not change with time 
and the equivalence class of the pole divisor of ip(Q) evolves linearly on the Jacobian of T. 
The algebraic-geometric procedure allows us to construct ip{Q) an d L(z) explicitly in terms 
of Riemann theta-functions, given an appropriate Riemann surface and a divisor on it. This 
procedure does not require any Hamiltonian theory, although, in most cases, the correspond- 
ing physical systems are governed by completely integrable Hamiltonian equations. 

Krichever and Phong [5j [7j suggested a general approach to the Hamiltonian theory of 
integrable systems with Lax-type equations. They introduced a two-form on the space C 
representing a Hamiltonian structure of the system. Their formula is universal and works 
even in the infinite-dimensional case. In our case, it is defined as: 

Un = ~\ Yl res Tr {^~ lLl ~ n5L A &® - A K l ~ n 5K) dz = — res nd z, (3) 

where \& is an eigen-matrix of L(z), i.e., L^f = and its columns are just vectors ip(Q) 
on different sheets of T. The sum is taken over the poles of L(z), the zeroes of det L(z), 
and the poles of dz (if there are any). The number n is an integer parameter. Formula 
is well-defined on the space £, but it depends on the normalization of if>(Q). It becomes 
independent of the normalization when restricted to the leaves, where the one-form k l ~ n 5kdz 
is holomorphic. As a by-product, ui n becomes non-degenerate and independent on gauge 
transformations L —>■ gLg~ l on the leaves (for the proofs see [9] and [10]). 

An alternative approach to the Hamiltonian theory of integrable systems uses a so-called 
r-matrix. An r-matrix defines Poisson brackets on £ between L(u) and L(v) for fixed u and 
v. In particular, when the domain of L(z) is a Riemann sphere, then, the simplest non-trivial 
r-matrix is 



V 



z 

t,3= 



We denote the domain of L(z) by V , and call z a spectral parameter. When L(z) is a 2 x 2 
(so N = 2) matrix function, an elliptic r-matrix (i.e., T is an elliptic curve) is equal (up to 
a scalar factor) to ([I]): 

1 B' n (6 01 {z) 9 00 (z) 9 w (z) \ 

2n U n (z) \ U 01 U 00 V 10 ) 



where cx 7 - are the Pauli matrices 



1\ fO -A I 



The notation for the Jacobi theta-functions Oij(z) is the same as in Mumford |14j . Notice, 
that in the 2x2 case (here do is the identity matrix): 

^ = \ ( a ® °"0 + °\ ® <*\ + 0"2 ® &2 + 0"3 ® ^3) • 
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These r-matrices satisfy the classical Yang-Baxter equation 

[r 12 (u - v), r 13 (u)] + [r 12 (u - v), r 23 (v)] + [r 13 (u), r 23 (v)] = 0, 



which holds in the space C <8>C <8>C with the evident notation T\ 2 = r®/, r 2 % = I®r, 
etc. 

For each r-matrix, one can construct two types of Poisson brackets: a linear bracket and 
a quadratic bracket. A linear (or Lie algebraic) bracket is defined as (PQ) 



These Poisson structures and the two-forms u n are degenerate on the whole phase space C. 
However, they naturally foliate the space into symplectic manifolds, where it is possible to 
compare them. 

The classical Yang-Baxter equation is a sufficient condition to ensure that the brackets 
above are indeed Poisson. There is no universal way to construct solutions to this equation, 
although many solutions are known (|15j). 

Ui, given by ([3]), coincides with the linear brackets. It is possible to show (|9]|) that the 
two-form u\ coincides with a Kirillov-Kostant form on the direct product of coadjoint orbits 
of the GL(N) action. When the poles of L(z) are simple, the corresponding principal parts 
can be identified with Lie algebra duals qI*(N). Symplectic leaves are determined by the 
condition that the one-form SKdz is holomorphic, or, equivalently, that the principal parts 
of K at each pole of K are fixed. The latter condition fixes some orbits in qV(N). 

It is shown in pU] that quadratic brackets coincide with u 2 . The nature of quadratic 
brackets is more complicated. Quadratic brackets in the elliptic case were explicitly com- 
puted only recently in [I7J for a multi-pole Lax function. Formulas (3.2)-(3.6) in [17] appear 
complicated with no apparent pattern. In this paper we show that if one uses a multiplicative 
representation for a multi-pole Lax function, then the quadratic bracket assumes a remark- 
ably simple form. In the case considered in [T7J , the bracket may be obtained as a reduction 
from the direct product of spaces with single-pole quadratic brackets. 

The simplest illustration of these ideas is a rational case, i.e., when T is a Riemann sphere 
and L(z) is just a meromorphic matrix function with a fixed number of poles on the extended 
complex plane. We assume that L(z) is in "general position," which means that it belongs 
to a big open cell in the space of meromorphic functions. In particular, L(z) may only have 
simple poles with residues of rank one, and it is diagonalizable at least at one point. Without 
loss of generality, we assume that it is diagonalizable at z = oo. Other cases, when L(z) has 
a higher-order pole or a residue of a higher rank, are special and may be obtained as a result 
of some limiting procedure. 

The most natural way to write a meromorphic function with simple poles is to specify 
positions of its poles and their corresponding residues. This leads to the formula 



{L(u) f L(v)}i = [r(u - v), L(u) <g> I + 1 ® L(v)), 



(5) 



and a quadratic (or Lie group) bracket is 




(6) 
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where a«, foj are iV-dimensional vectors and L is a constant matrix. We call it an "addi- 
tive representation," because it reflects the additive Lie algebraic structure on the space of 
Lax functions C. This representation is well-suited to the linear Poisson bracket (or the 
corresponding symplectic form u>i), which, in these coordinates, equals 

d 

u>i = 5af A 5bi. 
i=i 

However, different coordinates are natural for the quadratic bracket. Any function L(z) has 
an equivalent "multiplicative representation:" 

L = L (i+P^Jl+M-)^ 

V z-zxj \ z-z 2 J \ z-z d J 

where Pi and are also iV-dimensional vectors. 

It seems that the coordinates Pi, qi appear first in [6] , and later they are used by Borodin [12] 
in the theory of difference equations. These coordinates emphasize the multiplicative (or 
"group" ) structure on the space of Lax operators. Essentially, a multiplicative representa- 
tion is a particular case of an "integrable chain" from [10] . However, in [10] integrable chains 
were not related to an additive representation, and it seems that the name "multiplicative 
representation" is more appropriate than "chain" in the present situation. The quadratic 
2-form corresponding to the quadratic Poisson brackets equals 

d 

to 2 = ^2SpJ A6 qi . (7) 

i=i 

Notice that it is a highly non-trivial task to arrive at explicit commutation relations between 
coordinates Pi,qt (or, tjj and hi in the linear case) starting from Formulas (JHJ) or (JHJ) ■ One 
advantage of Krichever- Phong's Formula (j3]) is that it allows us to find them and that it 
works equally well in rational and elliptic cases. A multiplicative representation also exists 
in the elliptic case, and similar statements related to the quadratic bracket hold (see Sections 
4-6). We consider separately the case of "general position," when the poles of L(z) are simple 
and have rank one and the Sklyanin case of higher rank poles. 

It turns out, that for Sklyanin's Lax Matrix [2j, the 2-form u 2 given by Formula (E]) 
is degenerate even on the leaves where 5 In kdz is holomorphic. In order to circumvent 
this difficulty, we need to introduce an additional parameter a to 4 Sklyanin's variables 
so, Si, s 2 , S3. Surprisingly, this provides an extension of a quadratic Poisson algebra: 

{s , si} 2 = -8^s 2 s 3 , {s , s 2 } 2 = 6>o sis 3 , {s , s 3 } 2 = -6 ) i sis 2 , 

{Sl, S 2 } 2 = —SqS 3 , {Si, £3)2 = SqS 2 , {s 2 , 53)2 = — SqSi, 

{u,s } 2 = s , {u, s 1 } 2 = s 1 , {u, s 2 } 2 = s 2 , {u, s 3 } 2 = S3. 

These brackets coincide with the Sklyanin Brackets [2] without the last 4 identities. However, 
the symplectic leaves for the original Sklyanin brackets have dimension 2, whereas for our 
extension the leaves have dimension 4. The leaves are determined by one condition 

5((sl + sX + sjO^/isl + sj + si)) = 0. 
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As we said earlier, Formula ([3]) provides an hierarchy of 2-forms on the space of Lax 
operators C. Different integer values of n correspond to different 2-forms. In particular, it 
provides a non-trivial cubic bracket on the space of Sklyanin's Lax functions. The cubic 
bracket corresponds to u 3 and is given by relations 

{so, si} 3 = -26 , o 1 sos 2 s 3 , {s , s 2 } 3 = 29^ s s 1 s 3 , {s , s 3 } 3 = -2^ s sis 2 , 

{Sl, S 2 } 3 = S 3(-Si^00 + S 2^01 _ S o)> { S l> S 3}3 = S 2 (sl - s\0 1Q + S 3 Q1 ), 

{s 2 , s 3 } 3 = si(s 2 + s 2 2 9* + s 2 3 0% ). 

This bracket is non-degenerate on the leaves 5u = 5((sq + s\6qq + S2^oi)/( s i + s 2 + s§)) = 
S((sl + s| + s§)/so) = of dimension 2. 

The cubic bracket is related to the quadratic one in the following way: 

{L{u) f L{v)} 3 = ~[{L{u) f L(v)} 2 , L(u) ® <r + a ® L(v) - s ^o <g> <t ] + , (8) 
where [A, + = AB+BA. Notice, that the quadratic bracket also has the similar expression: 

{L{u) f L(w)} 2 = -[{L(«) f L(v)}i, L{u) ® a + a ® - s (Xo ® Oo]+. (9) 

It would be interesting to know the following: whether the last 2 formulas generalize to 
an arbitrary Lax function, whether the whole hierarchy of symplectic structures may be 
obtained in this way, and what implications it has on the integrability of Lax equations. 

2 Linear form in the rational case 

The coordinate form of Formula ((Sj) is 

{Lij(u), Li s (v)} = - j-^ ((Lij(u) - Lij(v))Sis + (L is (v) - L is (u))8ij) . (10) 
It is instructive to see that in the additive representation of the Lax operator 

L(z) = L + J2^-, (11) 

Krichever-Phong's Universal Form ([3]) corresponding to n — 1 equals 

d 

u 1 =^2SaJ A8b h (12) 
i=i 

and its inverse 

{^03 = ^ = i b l b j} = ( 13 ) 

agrees with Brackets ffTUl) . 
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Let us introduce matrices K and K\ as coefficients of the Laurent expansion of K(z) at 
z = oo 

K = K^ + o(^). 
We assume that L is in a diagonal form, i.e., L = K . 

Theorem 1. Universal Form corresponding to n = 1 and Lax Operator 07]) Ct>i = 
X)i=i A ^ S a well-defined symplectic form (independent of the normalization of \& 

and o/ gauge transformations) on the leaves 5Kq = 5K\ = and 5(bjai) = 0. 

Proof. In general, Form ([3]) depends on the normalization of First, we determine condi- 
tions that make u n independent of the normalization, and then we compute U\. 

A change of the normalization corresponds to a transformation \& — > tyV, where V is a 
diagonal matrix. Formula transforms as follows: 

-> £1 + 2Tr (K^dK A . 

If the one-form K x ~ n 5Kdz is holomorphic on T, then the second term in the last formula 
does not contribute to u n . Or, equivalently, we should restrict the one-form k x ~ n 5kdz to 
some leaves where it is holomorphic. 

Recall that n = 1 under the assumptions of the theorem. The one-form 5kdz may have 
poles on T above z = oo or above the points z = Zi (i.e., at the poles of k). Let ki/(z — Zi) be 
the principal part of k(z) at one of its poles and ipi{z) be the corresponding eigenvector of 
L(z). Since the principal parts of both sides of the equation Lip = kip must be equal to each 
other, we have aibjipi(zi) = ipi(Zi)ki. Multiplying both sides of the last equation by bj on the 
left and dividing by bfipi(zi), we obtain that ki = bjai. The one-form 5kdz is holomorphic 
at z = Zi when its principal part vanishes, i.e., S(bfai) = 0. Since dz has a second-order pole 
at z = oo, the one- form Skdz is holomorphic at oo if and only if 5K = 5K\ = 0. 

As a by-product, it turns out that u>i is symplectic on these leaves and does not depend 
on gauge transformations L — > gLg' 1 , where g G GL(N). See [H] for the proof when n = 1. 

Now, we evaluate Q. The second term in Formula (j3J) has poles only at z = Zi, z = oo, 
branch points of T, and at the poles of ^(z). When the eigenvectors of L(z) are normalized 
so that the sum of their components equals one, then \&(oo) is the identity matrix and 
<5\I / (oo) = SK (oo) = 0. One can check that the residues of Qdz at z = oo and at the branch 
points vanish. Since the sum of all residues of a meromorphic differential must vanish, we 
can rewrite as 

1 d 

ux = — VresTr (V^SL A dz. (14) 

2 z i 

i=l 

Clearly, 

resTr (ty~ 1 5L A 5V) dz = Tr (^- 1 (z i )5(a i bJ ) A • (15) 

Only one entry of the matrix function K(z) has a pole at z — Zi. Without loss of generality, 
we assume that its principal part is diag(aj, 0, 0)/(z — zi). The identity = ^/K implies 
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that 

aibj^(zi) = #te)diag(oi,0,...,0) and *- 1 (^ i )o i 6f = diag(a i , 0, O)*- 1 ^). 
Consequently, 

6f = (A, 0, 0) and ^(z^ = (A, 0, .., 0) T . 
From the last two identities we deduce that bfdi = (3if3i and 

bJd^iz^-^Zi) = bj5\n(3i-5bj, S^iz^izi)^ = 5a { -a^ln/3,. 

Therefore, we deduce that (1T51) equals 

resTr (* _1 (JL A tfc = 256f A (Jo*, 

which completes the proof. 

□ 

We now check directly that Brackets (TT5|) agree with r- matrix Poisson Brackets (jlUp . 
Using the properties of tensor products, one can show that 

{a t bj f a^f} = (a; ® /) (/ ® bj) -(I® a t )(bf ® /) = (a^f ® I)P - ( J <g> a^f )P. (16) 

Clearly, {a^fef f cijbj} = for % ^ j. Formula (J5J) follows if we consider L as a constant 
matrix and use the properties of the permutation matrix B ® A — V(A ® I?)? 7 and V 2 = I. 

A Lie algebraic interpretation of the linear brackets has been suggested in [9]. For- 
mula (I14p may be rewritten as 

d d 

Ul = V res Tr (LcW^ 1 A cW^ 1 ) tfe = V c^. 

8=1 m=l 

Let us define L m as 

L ( z ) = + 0(1) = + 0(1). 

2r Z m Z Z, m 

Then we can identify L m with a point of gl*(N) and the Lie algebra with its dual using 
the Killing form. Each term uj' m equals the Kostant-Kirillov form defined on an orbit of a 
co-adjoint representation of a Lie group. As we saw before, 5(6^a m ) = on the symplectic 
leaves, which corresponds to the choice of some orbit in the Lie algebra. Therefore, u)\ is the 
Kirillov-Kostant form on the direct product of d coadjoint orbits of GL(N). 

The Poisson brackets that correspond to each Kirillov-Kostant form are 
and they coincide with Formula ([161) . 
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3 Quadratic form in the rational case 



As stated in the introduction, a rational matrix function L(z) in general position with d 
poles has 2 equivalent representations: 

ab T 

• an additive representation L(z) = L + J2i=i — ~~ L ~ > anc ^ 

Z Zj 

• a multiplicative representation 

L = lJi + M.) ( I + ™L) Jl+ ™L) = W,..B„ 

V z-zij \ z-z 2 J \ z-z d J 

where Oj, fej, p,, and are iV-dimensional vectors and Lq is a constant matrix. 

On the symplectic leaves, Formula ([3D is invariant with respect to gauge transformations 
L — > g~ 1 Lg, — > (7 1 \I/ , where g e GL(N). Therefore, we may assume that the matrix Lo 
is diagonal. 

The following lemma proves the equivalency of additive and multiplicative representations 
for arbitrary d. A similar result has been proved by Borodin in |12j . Dzhamay [18] has proved 
the equivalence when d = 2. 

Lemma 1. For any meromorphic matrix function L(z) corresponding to a Zariski open 
subset of parameters (zj,aj,6j), there exists a multiplicative representation. The converse 
also holds. 

Proof. An additive representation follows immediately from the multiplicative one by taking 
the residues at the points Zj. 

To prove the converse, we assume that we have an additive representation and construct 
vectors p, t and gj. Let z± , z%, z~[ be zeroes of det L. Notice, for a multiplicative repre- 
sentation, one has z~ = Zi — pjqi- Let ip* be a left eigenvector of L and the corresponding 
eigenvalue k have a pole at Zd- If the principal part of k is C/(z — zj), then the principal 
parts of both sides of the equation ip*L = kip* are 

^*{z d )L B 1 (z d )...B d . 1 (z d )p d -^— = if>*{z d ). 

z- z d z- z d 

Since ip*(z d )LoBi(z d )...B d ^i(z d )p d is a number, the latter equation implies that ip*(z d ) oc gj. 
Likewise, if ip is the right eigenvector L~ Xr p = k~ l ip, where k has a zero at z = z J", then 
ip(z^) oc p d . Since PdQd = Zd — z^ , we can recover pd and qd up to a scaling factor, which 
does not affect Bd- We can repeat this procedure for the conjugated matrix BdLB^ 1 to find 
B d -\. We can find all factors B\, B2, ...Bd in this manner. 

The only thing left to prove is that BiB 2 ...BdL~ 1 is a constant matrix. By construction, 
L(z d ~)pd = 0. By assumption, L(z) is in general position, which means that dimker L(z~[) = 1 
and that the residue res^- L~ l (z)dz has rank 1. Since LL~ l = I, it must be that 

L -i {z) = i^L + 0[l) 

Z ~ Z cL 
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for some vector v d . Since B d (z d ~)p d = 0, we conclude that the function B d L~ x is holomorphic 
at z^. By construction of the vector q d , we have q r ^L~ 1 {z d ) = 0, which means that B d L~ x is 
also holomorphic at z d . Using the same arguments, we can show that 

B d L-\z) = ^Ll + (l) 
z ~ z d-i 

and that q^^B^Zd-ijL" 1 (z d -i) = 0, which implies that B d _\B d L~ x is holomorphic at z~^_ x 
and Zd-i- By induction, we prove that BiB 2 ...B d L~ l is an entire function on the extended 
complex plane, hence it has to be a constant. □ 



Now, we are in a position to prove: 

Theorem 2. Universal Form (Ejj corresponding to n = 2 and a rational Lax matrix L(z) = 
L BiB 2 ...B d in the multiplicative representation equals 

d 

uj 2 = ^5pJ Adqi. 
i=i 

Symplectic leaves are determined by the conditions SKq = SK\ = and S(qfpi) = 0, where 
K(z) = Ko + Kt/z + Oil/z 2 ). 

Proof. The proof is similar to the proof of Theorem [TJ We can rewrite Formula ([3]) as 

1 d 

u 2 = — V res TV (V^L^SL A SV) dz, (17) 

where z^ = Zi — qfpi are zeroes of det L. 

The two-form uj 2 is symplectic, independent of the normalization of \l/ and of gauge 
transformations, provided that the one-form S In Kdz is holomorphic on V or S In kdz is 
holomorphic on V. Since dz has a second-order pole at z = 00, we should fix Kq and K\, 
i.e., two conditions that determine symplectic leaves for uj 2 are SKq = SK\ = 0. Other 
possible singularities of Sin kdz are at the points Zi and z~ . One can check that Sin kdz is 
holomorphic if Sz; t = Sz^ = 0, which yields the condition S(qfpi) = 0. 

Let us introduce matrices T d = L B 1 B 2 ...B d , T d _ x = B d L BiB 2 ...B d ^ 1 , 7\ = B 2 B 3 ...B d L 
Since T d = L, we have T d ^/ d = ^ d K and = Matrices Tj with i < d are conjugated to 
T d , i.e., = *tK, where ^ = ^ d _ 2 = B d ^B d ^ d , * x = B 2 B 3 ...B d V d . 

The following transformation of co>2 is almost identical to the one used in [TU|- One can 
show that: 

d 

Tr {^! d l T d l ST d A 5tf d ) = ^ Tr (^fB^...B^SB k B k+1 ...B d A 5tf d ) = 

fc=i 

= ]TTr (^S^* A " E Tr {B^-B^SBm A *(B w ..-5d)) ■ 

fc=i fc=i 
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Notice that the last sum does not have any poles except at the points and z i and vanishes 
after the summation over all residues. Therefore, 

1 d 

U 2 = --J2 reS Tr (*J lB j l5B 3 A ^j) dz - 



The matrix ^ consists of normalized eigenvectors and it does not have poles at the points 
Zi, z~ for any % in general position. However, matrices tyj may acquire poles at the points 
Zi,z~ for % > j. Since matrices tyj (j < d) consist of eigenvectors of Tj, we can normalize 
them: tyj = ^jFj. The matrix functions F^ 1 are diagonal, possibly having poles at Zi,z^ 
for i > j. Normalized matrices ^ are holomorphic at z^, z~ for any i. 

The second term on the right hand side of the identity 
Tr '/; ; dllj A 8Vj) = Tr U>^ 1 Bj 1 8B j A 5%) - Tr (v^BjHBjVj A 6\nF^ 

is holomorphic at z iy z~ for i > j, because 5zi = 5z^ = 0. 
Therefore, our formula for u 2 becomes: 

1 d 

u 2 = — res Tr {^B^SBi A 5%) dz. 

2 i=l **'*< 



Plugging in the expression for Bi and computing the residues, we obtain: 



u 2 



d 

-5E 



QiPt 



AWi(z~ 
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Let us fix an integer i. Define the function Uj as Tj = L^Sj. The function C/j is holo- 
morphic at Zi, and Bi has a simple pole there. \I/j is holomorphic at Zi, and i-T is a diagonal 
matrix with all but one entry holomorphic at Zj. Without loss of generality, assume that 
Kn has a simple pole at z\. The principal part of the identity UiB^i = ^>iK implies that 
Ui(zi)piqf^i(zi)^~ (zi) is a diagonal matrix and qf^i(zi) = (««, 0, 0, 0), where a« is some 
scalar function. Taking the variation of the latter identity, we deduce: 

qj5^ l (zi)^; l (z i ) = qJSlnai - 5qf. (19) 

Similar arguments for ty^Tf = K^ 1 ^^ 1 at the point z^ prove that \E^~ {z^)p% = (Pi, 0, 0) T 
and 

8%{zr)^7\ z -) Pi = Spi-piShiPi. (20) 

Substite f[T9~|) and f[2"Uj) into (|18p to complete the proof of the theorem. □ 
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We now check that uj 2 from Theorem [2] agrees with r-matrix Brackets (Q. The Poisson 
brackets corresponding to u 2 are 

{Vii ? rf} = (Pi ® /)(/ ® <?f ) - (/ ® Pi)(gf ® /)■ 

One can check that 

{B , ( u) ? *(„)} = fe^(£MV(£^MM. 

Consequently, we deduce: 

{Bi(u) f = <%[r(u - «), <g> ^-(t;)]. 

Then, if we consider L as a constant matrix, the group property of the quadratic bracket [16] 
gives Formula (jSJ): 

{£(«) ® L(v)} = [r(« - v), L{u) ® L{v)\. 

4 Elliptic case: general position 

Certain difficulties arise if one wants to construct a non-trivial Lax Equation ([I]) on an elliptic 
curve. In this case, the principal parts of an elliptic function L(z) can't be arbitrary due to 
the relation 

res L{z)dz = 0, 

i 

which is not invariant under Flow ([I]). 

Two general approaches are known to overcome this difficulty: one of them is due to 
Krichever and Novikov [3], another one is implied in Sklyanin's Paper [2]. We consider 
Sklyanin's approach in Section O The idea of \4\ is to introduce N additional poles to the 
functions L(z) and P(z) with special dependence on t, so that Equation (JTJ is non-trivial. 
Positions of the "main" poles and their principal parts, i.e., the set of data (z iy res 2i L(z)dz), 
determine the function L(z) up to a complex scalar. In order to avoid "pathological" cases 
(e.g., when some of the poles Zj coincide), we consider the values only in a Zariski open 
subset. We denote the divisor of the "main" poles by D + . In the same way as in the rational 
case, we assume that all poles are simple and the residues at the points 2, have rank one. 

The additional poles with coordinates q\, 1 < % < N play the role of so-called Tyurin 
parameters that parametrize framed stable degree iV holomorphic bundles. We impose the 
constraints: 



L(z) = \ + Ln + o(l), ajfii = 0, ajLn = aj^ 
1i 



where the matrix (cti, a 2 , a^) is the identity matrix. If V 1 is a bundle that corresponds to 
(q{, then L(z) may be identified with a section of H°(T, End(V l )(D + )), where the order 
of poles is bounded by the divisor D + . Lax equation defines an evolution of the bundle 
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V 1 . This construction is also applicable to Lax equations on a Riemann surface of arbitrary 
genus (see [H] for details). 

The linear symplectic form u\ was computed in |9J. In this section we compute a quadratic 
form for a multiplicative representation of L(z). It turns out that, in order to define a 
multiplicative representation, we need to introduce a sequence of vector bundles. Suppose 
that we have a representation: 

L(z) = B d B d ^...B x . (21) 

The data (zj, res^ L(z)dz) depend on 2dN parameters. In order to completely determine the 
function L(z), we must specify ^2^ =1 q[- Therefore, L(z) depends on 2dN + 1 parameters. 
If we assume that each function Bi is holomorphic at Zj for j ^ i, then the right hand 
side of (T2~T1) depends on (2iV + l)d parameters. However, the problem is that the Tyurin 
parameters for each Bi are different, i.e., functions Bi are endomorphisms of different vector 
bundles and their composition does not make sense. 

A way around this difficulty was suggested in [TD] . The idea is to consider a sequence of 
vector bundles V m corresponding to the Tyurin parameters (q l m ,ai), such that V d+1 = V 1 . 
Then each B m e Hom(V m , V m+1 )(z m ) and has iV poles at the points q l m , so that 

Bm = ^K + 0(1)t (22) 

Z 1m 

where the f3i(m) are iV-dimensional vectors. The inverse functions are homomorphisms of 
vector bundles, and the vector bundles are in the opposite order: 

where z~ is the pole of L' 1 ^) which is distinct from the points q l m+ i and D_ = ^2 
The inverse functions also have poles at the points q m+ i, 1 < i < N, so that 



m m 



fi _ 1 = « +0(l) (23) 

The function B m (z) is elliptic with a simple pole of rank one at the point z = z m and may 
be written out explicitly using Weierstrass sigma functions: 

m a(z- z m )a(q l m+l - q m )a(z - q> m ) 
Its inverse has a simple pole of rank one at the point 

N 

Z m = Z m + ^ ] {q m — Q m +l) ' (^^) 
i=l 

where the complex numbers q l m are subject to the periodicity conditions q l d+1 = q{. 

Since the transformations f l m — > X m f m don't change L(z) provided that rim=i ^™ = 1, 
the total number of independent parameters needed to describe a chain Bi, B 2 , B d is 
2dN + 1, which coincides with the number of parameters for the function L(z). 
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Theorem 3. The quadratic symplectic form (given by ^)) equals 



^ = 1^ 1^ 6 ln FT -a p ^ I A 6<lm - (25) 

m =li=l \ Hp=l°Wm % 



m+1 / 

Symplectic leaves are determined by the conditions 5z m = 5z^ = 0. Only 2d — 1 of them are 
independent and in general position, so the dimension of the leaves is 2d(N — 1) + 2. 

Proof. As opposed to the rational case, Qdz in Formula ([3]) has poles at the points q\. 
Moreover, the residues of the two terms of Qdz are not equal to each other. We can rewrite © 
as: 



1 1 - 

u 2 = — res Q'dz + - ^ res A K'HK) dz, (26) 

2 2 i=i 95 

where f2' = Tr (ty~ 1 L~ 1 5L A and the first sum is taken over the points q\,z m ,z^, 1 < 
i < N, 1 < m < d. 

Consider the following matrices: T\ = L = B d B d _ 1 ...B 1 ,T 2 = B 1 B d B d _ 1 ...B 2 , ...,T d = 
B d -\B d _2---B 2 B\B d . Let \I> m be an eigen-matrix of the corresponding T m , i.e., T m ^ m = 
ty m K. Matrices with z > 1 are related to = \& in the following way: ^ 2 — -Bi^i, ^3 = 
S a Bi* 1 ,...,* d = J B d _iS d _ 2 ...Bi* 1 . 

One can rewrite fi' as: 

d d 

fi' = ^ Tr {^B^8B m A 5tf m ) - ^ Tr (B^B 2 l ...B^8B m A <f(B TO _ 1 B m _ 2 ....B 1 )) . (27) 

m=l m=2 

In general position, fi'cfe does not have any residues at the points q l m with m > 1. Therefore, 
we can safely add them in the summation in Formula (126]) . However, the second term in (1271) 
vanishes after this procedure. The first term in Formula (126]) becomes: 

1 1 d / JV TV \ 

~2 ^2 Tes ^' dz = ~2 ( res ^ fi ™^ + r es 2 -^ m ^ + ^res 9 j n fi m rfz + ^ res^ Q m cfe J , 

m=l \ i=l i=l / 

(28) 

where Q m = Tr {m^B m 1 5B m A 5* m ) . 

The functions \l/ m may have some poles at the points and with i < m, but an argument 
identical to the one used in the proof of Theorem [2] shows that they don't contribute to uj 2 . 
Now we compute each term in Formula (j2"8|) : 

One can check that B~} is (the computation of L" 1 in [8J is almost identical to this 
computation): 

l B -\\ = °( Z ~ Z m + <li- gm+lM* ~ 1m) x Up=l ^(iL+l ~ 1m) Up^l °~(gm ~ C+l) 
M fL<T( z ~ z m)°( z ~ lln+l) Ylp^k a (lm ~ Qm) Up^l ^(iL+l ~ Im+l) ' 
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Below we will use the following notation: res Zrn B m dz = PmQ^, res z - B^dz = P m Q m , where 
P m , Qm, Pm-, Qm are the ^-dimensional vectors: 



FT a(n k - n p ) 
Pa(z -a 1 ) P k - cy(z- a k ) p=1 — 

Jm u \ 4, rn H m +l)i 1 m u \^m Hm) 77 „(„k _ rf \ 

Jm) 



, , j s p^i 1 rip=l ^(Qrn+l 9m) 



fm a ( Z m 9m+l) Ylp^l a ( ( lm+l 9m+l) 

We now compute res Zm Q m dz. Since T m has a simple pole of rank one at z m , one diagonal 
entry of K has a simple pole at that point and \l/ m is holomorphic at z m . The principal part 
of the equation T m ^/ m = ^ m K implies that Q^$ m (z m ) = (x m , 0, 0). Taking the variation 
of the last identity, we deduce that: 

Q^S^mizm)^^ 1 ^) = Q T m b\wH m - bQ T m . 



-l 

m 



Using the identities i? m 1 (2; m )P m = and Q^B^^Zm) = (which follow from B m B, 
B^B m = I), one can show that: 

res Zm tt m dz = bQ T m A B m 1 (z m )5P m . 

The identity ^^T~ l = at the point £~ implies that \E f ~ 1 (z~)P m = (n m , 0, 0) T 

and, consequently: 

5^ m (^ m )^ m 1 (2; m )P m = 5P m - P m 5\\iK m . 
Using the fact that Q^Pm^m) = B m (z m )P m = 0, one can prove that: 

ies z -Q m dz = -5Qm A B m (z m )5P m . 

Now, we compute the term res^fim*^ in (J2SJ). We assume that the poles of \l/ m do not 
coincide with any of the points q l m . Then the Taylor expansion of is: 

m m = tf m (<4) + *' m (Qm)(* - <L) + ^m(?m)(* ~ 9 m ) 2 + 0((* - gj 3 ). 

The principal part of T m ^f m = ^ m K implies that aj^ m (q l m ) = since ^ m K is holomorphic 
at q m . However, T m = ^/ m K^^ has a pole at q m , so the Laurent expansion of ^P" 1 is: 

*i? = + < + *m(^ - 9 j + 0((* - 9 m ) 2 )- 

Plugging the series of \I/ m and ^l/" 1 into ^m^" 1 = J, we obtain: 

^ m {q m )li{m)aJ = 0, * m (^)* m + *' m (q m Mm)aJ = I, (29) 

Formulas (j22"]l. (|25|) and the fact that B m B^ = B m 1 B m = I imply that: 

ajB m (q m+1 )=0, ajB-\qi n )=0, B m \q j/3;(m) = 0. (30) 
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Now, we show that: 

Clearly, we have: 
Using (1221) . we obtain: 



Q; 



dz 



AM = 1. 



(31) 



dz 



dz 



AM 



AM 



'I'm 



-«jB-^B-%(m). 



d{B^] 



dz 



AM 



so (HP follows. Now, using Formulas p2]h (|23]|, (J29])-([3ID, we compute: 

r ^( 5 m' 



res„ 



(JAM A <5<4 + a/ ^(^^(m) A &4. (32) 



In order to compute res g ; Q m dz, we need 2 additional identities that follow from \l/ m+ i = 
B vl/ ■ 

*m(9m+l)AM = 7i(™ + 1), ^mW m+1 H(m + 1) = A(m). (33) 
Using ([23]), (J29])-([30D,([33D, we show that: 



res„i Vtmdz = a, 



dB rt 



dz 



8Pi(m) I A 5ql n+1 - aj ^' m+1 {<L + x)H{m + 1) A 5q l m+l . (34) 



Notice that the last terms in ( 1321) and ( 13~4"1) telescope after summing over i. The only terms 
that we need to be careful about are those with i = d + 1, since ^d+i = ^ One 
can check that aj%(q\)6^i(l) = af^' d+1 (q d+l )5-fi(d + 1) + af% (q[)5 In K(qi)%(l) and 

res (V^SVi A K- X 8K) dz = af^f' 1 (qi)5]iiK(qi)^ i (l) A 8q[, 

so we finally arrive at the formula: 

d 



"a = - ~ £ [ 5 Q™ A B m \z m )5P m - 5Ql A B m (z-)5P m + 



(35) 



m=l 



N 



dB n 



<9z 



A' 



<?m+l 



<5A(m) A <^ +1 - 



a 



i=i 



dz 



5AM 



Straightforward, but rather lengthy computations show that: 



d N 



uj 2 



fL-i U P ^i °{<L - <il 



p > 

m+li 



A5q^ 



on the leaves, i.e., when tfz™ = 5z m = 0. 



□ 
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5 Sklyanin's case 



Sklyanin [2] has defined quadratic Poisson brackets on the space of Lax matrices with one 
pole: 

L = «/ + , *) + (i±I, r) + g, ,) , (36) 

where 



9ii(to|r)9i,(z|T) 

and #ii(z|t) is a Jacobi theta function. The notation for Oij(z) is the same as in Mumford 
The functions 8ij(z) and 9ij(z\r) correspond to the curve C/(Z+rZ). The function 9ij(z\2r) 
corresponds to C/(Z+ 2rZ). 

The function L(z) has the following translational properties: L(z + 1) = a^L{z)a 3 , and 
L(z +t) — <j\L(z)<J\, where (7j are the Pauli matrices: 

l\ fO -i\ (I 



a> \i o)^ 2= \z o J^ 3 vo -i 

The function L(z) is elliptic on the curve C/(2Z + 2rZ), and, due to the prescribed mon- 
odromy properties, the sum of its residues is automatically zero. This construction allows us 
to choose the principal part of L(z) at the point z = arbitrarily provided its trace vanishes. 

L(z) may be identified with an endomorphism of a vector bundle on the elliptic curve 
C/ (Z + rZ) with a pole at z = 0. This bundle has degree 1 and rank 2 and is described by 
its section s = (si(z), ^(z)), that transforms according to the formulas: 

s T (z+l) = Qs T (z), and s T (z+r) = As T (z), where Q = cr 3 , and A = exp (— iri(z — U\) — mr/2)ax. 

Here, we are following the notation from [IT] . One difference with [17] is that we have 
introduced an additional parameter u±, which changes the vector bundle, but does not affect 
L(z). We need it to obtain proper symplectic leaves for the 2-form ([3]). 

Formula ([3]) has been proved to work only in the case of degree 2 bundles, but we can 
still apply it to Sklyanin's case if we consider L(z) on the curve C/(Z + 2rZ) instead. As 
a result, the degree of the corresponding vector bundle doubles. A degree 2 bundle has 2 
holomorphic sections, and one can check that L(z) = gL(z)g^ 1 is an elliptic function on 
C/(Z + 2rZ), where 

'0oo(*-Ui|2t) 

9 w (z-u 1 \2r) / 

and u\ is an arbitrary parameter. Below we use letter u instead of u\ for simplicity. 

The function L(z) corresponds to a degenerate case of (!25|) . since its residue at z = has 
rank 2, whereas all residues in the previous section had rank 1. 

Theorem 4. Krichever-Phong's Universal Formula |3|) defines an hierarchy of symplectic 
structures on the space of matrix functions L(z) which, in turn, depend on 5 parameters 
(u, s , Si, s 2 , S3) . These structures vanish for n > 3, and uji,uj 2 , u 3 correspond to the following 
Poisson brackets: 
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• A linear bracket (n = 1): 

{si, S 2 }l = -S3, {si, s 3 }i = s 2 , {s 2 , S 3 }i = -Si, {«, S }l = 1, 

{w, Sj} = {s , Si} = for i = 1, 2, 3. 

27ms bracket is non- degenerate on the leaves 5u = 5s = 5(sl + s^ + S3) = 0, which 
have dimension 2. 

• A quadratic bracket {n = 2) : 

{■So, = —O01S2S3, {s , S 2 } 2 = Oq SiS 3 , {s , s 3 } 2 = — 6^ SiS 2 , 

{si, s 2 } 2 = -s s 3 , {si, s 3 } 2 = s s 2 , {s 2 , s 3 } 2 = -s si, (37) 

{U, S } 2 = S , {«, Si} 2 =Si, {w, s 2 } 2 = s 2 , 

{u, s 3 } 2 = s 3 . 

T/ie symplectic leaves for the quadratic bracket are defined by: 

s((sl + sle^ + sle 4 m )/(sl + sl + si)) = 

and have dimension 4- 

This bracket coincides with Sklyanin Brackets JE/ after the symplectic reduction to 
submanifolds u = 0. The submanifolds have dimension 2 and are given by the equations 
5u = 5(sl + s 2 2 + si) = 5(s 2 + s?£ 4 + sle 4 m ) = 0. 

• A cubic bracket (n = 3): 

{ s o, s i}3 — _ 26 ) q 1 s s 2 s 3 , {s , s 2 } 3 = 2#q s SiS 3 , 

{s , s 3 >3 = -20i O s o sis 2 , {si,s 2 } 3 = s 3 (s19q + s^i - «o), ( 38 ) 

{si, s 3 } 3 = S 2 (Sq - S^^o + S^gJ, {s 2 , S 3 } 3 = -Si(Sq + S^i + S^go). 

This bracket is non- degenerate on the leaves 

6u = 6((sl + s\9t Q + slel)/{s\ + sj + si)) = 5((sl + s\ + s^)/s ) = 
of dimension 2. 

Proof. Formula ([3]) is equivalent to: 

u n = k^mSk^) A Sztfi), (39) 

where ji are poles of eigenvectors ip of L(z) on f due to the following normalization: ipx = 1. 
In this case, the poles are given by the equation Li 2 (z) = 0. The symplectic leaves are 
determined by the condition that the one-form k x ~ n b~kdz is holomorphic. 

The proof of the theorem is a direct computation, and is similar in all three cases, i.e., 
when n — 1, 2, or 3. We outline it below in the case n = 2: 
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Points where the one-form k~ l 5kdz may fail to be holomorphic are zeroes and poles 
of k. They correspond to zeroes and poles of detL(z). The form will be holomorphic, 
if their positions are fixed, i.e., we need to impose the constraint 8zq = 0, where z is a 
zero of det L(z). When n > 3, uo n vanishes after we impose all necessary constraints. The 
determinant of L(z) equals: 



Hi 

therefore the symplectic leaves are determined by only one condition: 



det l(z) = 4 + s{et, + 4e* 01 + JH^M + s 2 2 + si), 



In the case n = 2, Formula (I3"9|) becomes: 

u 2 = J2 S]:ak (li)^W%)- (41) 

i 

The spectral curve is a 2-sheeted cover of the elliptic curve C/(Z + 2rZ) and if) has 4 poles. 
Two of them are located on both sheets above the point z = u+t+1/2. The other 2 are above 
the points z = z, z = 1 + 2t — z (where L 12 {z) = 0) and correspond to k = L 2 2(z) = L 22 (z). 
Sum (T4T1 equals: 

uj 2 = 5\nL 22 (z) A 5z - 5\nL 22 (l + 2r - z) A 5z + 5\ndetL(u + r + 1/2) A 5u. 
Equation (140]) implies that: 

5 In det L(u + r + 1/2) A 5u = 5 In (sj + s 2 2 + s 2 3 ) A 5u. 

Further computations show that: 

s 1 5s 2 A5s 3 -s 2 5s 1 A5s 3 + s 3 5s 1 A5s 2 2 2 2 

uj 2 = 2 -r-2 2 2\ \- 6\n(s^ + s 2 + s 3 ) A 6u. (42) 

s (s x + s 2 + s 3 J 

This form has rank 4 and corresponds to the Poisson brackets: 

{ s 0> s l} — — ^01 S 2"S3, {s , S 2 } = 6*o SiS3, {s , S 3 } = — ^ SiS 2 , 

{si, s 2 } = s s 3 , {si, s 3 } = s s 2 , {s 2 , s 3 } = -s si, 

{u, s } = s , {m, si} = si, {u, s 2 } = s 2 , {«, s 3 } = s 3 . 

Notice that the direct inversion of Formula (1421) leads to an additional factor 1/2 in all 
Poisson brackets, e.g., the first bracket is {so, si} = — 6 l o 1 s 2 S3/2. This factor appears because 
we double the elliptic curve, but we omit it in all formulas for Poisson brackets. 

Quadratic brackets for Sklyanin's case were also computed in [17] . and they coincide with 
those in [2]. Formula A. 23 in [17] yields the identities: 

* 2 (iW0=*X, eJI)- e JI^-)=«x, Ejl±r)- E2 (r^^ m , 



that allow us to simplify formulas in [T7] to: 

{s , st} = 27r0o 1 s 2 S3, {s , s 2 } = -2tt^ siS 3 , {s , s 3 } = 27r^ siS 2 , 
2 2 2 

{Sl, S2} = — {Si,S 3 } = S S 2 , {S2, S3} = — SqS].. 

7T 7T 7T 

The latter formulas agree with ours up to a constant factor after the rescaling so — > —sq/tt. 

Now, we have to compare the conditions which determine the symplectic leaves in [2] and 
in [T7] with fHUl) . The brackets in [2 J and in [17] have rank two, because there is no generator 
u. The symplectic reduction of our two- form (|42p to the submanifolds with the constant u 
yields 2 additional constraints (3 in total): 

5u = gtf + s l + = 6{s l + s 204 q + s 2 4 qi) = q 

The latter formulas coincide with the conditions for the symplectic leaves in [2] and in |17j . 

□ 

Remark: The Jacobi identity for Quadratic and Cubic Brackets (|37|) and (|38|) is equivalent 
to the only relation between 9 00 , O i, ^10; which is = ^ + 6*f . Therefore, one can get 
a 2-parameter family of quadratic and cubic Poisson algebras by replacing 9 01 and 6% by 
arbitrary complex numbers. 

The proofs for Formulas (IE]) and ([9]) are a direct computation using Riemann's theta 
relations. 



6 Degree 1 bundles with an arbitrary number of poles 

Sklyanin's brackets may be generalized to the case when L(z) has an arbitrary number of 
poles and any rank. An explicit computation was performed in |17j . In this section, we 
introduce a multiplicative representation for a multi-pole Lax function and show that it is 
natural for the quadratic brackets. For simplicity, we consider only rank 2 bundles. 

The construction of a vector bundle is the same as in the single pole case, and the Lax 
function L(z) with d poles has the form |17j : 



i=i 



s 1 / - 

Liz) =s°I + V &E x (z - Zi)I + ^ ai e m{z - z ^(f) ( - 
^— i J txi V2 

'1+T 



Zj) + 



s 2 

Til 



2 - Z -^) + S ^ 



- , Z Zj 



where the parameters s° satisfy Ylj=i s< j = 0; an d E\{z) = d z In 9u{z\t). 
The function L(z) has the following translational properties: 



(43) 



L(z + 1) = a 3 L(z)a 3 and L(z + r) = a\L{z)ai. 



(44) 
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It is an elliptic function on the curve C/(2Z + 2rZ), and, due to the prescribed monodromy 
properties the sum of its residues is automatically zero. This construction allows us to choose 
the principal parts of L(z) at the points z = Zj arbitrarily, provided that their traces vanish. 
The number of independent parameters needed to describe each function L(z) is 4d. In the 
same way as in the single pole case, the function L(z) becomes elliptic on T = C/(Z + 2rZ) 
after the conjugation L = gLg^ 1 , where 



9W 



6 m (z-u\2r) 

e w (z - u\2r) J ' 



The construction of a multiplicative representation is similar to the case of "general 
position." We introduce a sequence of degree 1, rank 2 vector bundles V m on C/(Z + rZ) 
described by their sections s(z), such that 

s T (z + 1) = a 3 s T (z) and s T (z + r) = exp (— m(z — u m ) — mr/2)aiS T (z). 

The factors B^ of a multiplicative representation are homomorphisms Hom(V\ V t+l )(zi) with 
a possible pole at the point z = Zi. We assume that V d+1 = V 1 . The following theorem 
relates additive and multiplicative representations: 

Theorem 5. Function L(z) has a multiplicative representation L(z) = B^Ba-i-.-Bi for 
a Zariski open subset of parameters P , where 



B m =s° m I<KA m , 4) + ^ie^<f> ( T - + A m , z' n 



+ S ^a 2 e^<j> + A m , A + ^a 3 <j> (\ + A m , z' m , . 

m \ 2 J hi \2 

A m = u m+ i — u m , z' m = z — z m , U\ = u d+ i, and the parameter u\ may be chosen arbitrarily. 

Proof. The function det L(z) is elliptic on the curve C/(Z + rZ) and in general position it 
has 2d distinct zeroes. Let z± , z 2 , be any d of them and we denote the rest with letters 
Zi,z 2 , ■■■jZj. Since det L(z) is an elliptic function, it must be that 

d d 

i=l i=l 

which makes possible to choose parameters u 2 ,u 3 , ...,«<* for any choice of U\, so that z~ = 
2z m — 2A m — 2~ for any m. 

Let us define the vectors P m , Q m , P m , Q m as: 

p o T P b T 

B- 1 = ±™% + 0(1), B- 1 = ±™% + 0(1). 



Then 



P O 



. -1 

(detB r 



in 



dz 



1 

-B 21 (z-)/B 22 (z-) 



(B 22 (z~) -B 12 (z~)) 
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PmQr, 



d 

dz 



(det B n 



~B m ( Z m )/B m (z n 



Therefore, we conclude that: 



\B m {z m ) B m (z m )*j 



( n 0ii 

rrst 



-Pm, 0^ 



+ A r 



V 



■is: 



#11 (A m ) 
i #01 [z m + A r 

#01 (A m ) 



+ IS 



9io(z m + A m ) \ 

m #lo(A m ) 

m / 

#00 (A m ) / 



and 



/ o 011 + A, 



Prn OC 



-7TS 



" #11 (Am) 
#01 [Zm + A 

Til J 



+ ^s 



#oi ( A r 



3 #10 (^m + A m ) \ 

m #io(A m ) 

#00 {Zm + A 

#00 (A m ) / 



where z„ 



>d-i- 



(45) 



..Bi, then Pi and Pi 



Notice that if one already has a representation = P^Pd 

span kernels of L(z) at the points 2 = z± and z = z[ . This is equivalent to two equations 
in si, s\, sf, sf. Using an additional constraint det Pi (z^) = 0, we can recover the values 
up to a common scalar factor for any function L(z) ([4"3]h We apply the same 



procedure to the conjugated functions BiLB 1 , B 2 BiLB l l B 2 , ... to construct d matrix 
functions Pi, P 2 , B d . 



Now, let us show that LB^ 1 B^ 1 ...B^ 1 is a constant matrix proportional to /. Since 
L(zi)P\ = 0, the function LPf 1 is holomorphic at z — . Likewise, it is also holomorphic 
at z — z^ . By construction, P ] ^ 1 (z 1 ) = 0, so LP-f 1 is also holomorphic at Z\. In general 
position, Pi (z 2 ) is a non-degenerate matrix. Therefore, by the construction of P2, it must be 
that L(z 2 )B^ 1 (z 2 )P 2 = 0, and, consequently, LB^ 1 B 2 1 is holomorphic at z± and z 2 . Using 
similar arguments and the fact that P~ 1 (^m) = for all m, we can show that LB^B^ 1 ...Pj 1 



is holomorphic everywhere on C/ (Z + rZ). 

Each matrix function B m is elliptic on T 
g m+ iB m g^-, where 



C/(Z + 2rZ) after the conjugation B n 



9m{Z) 



0oo(z 



u, 



J2r) 







# 



iolz 



|2r 



Using this fact, one can check that the product B^Bd-i-.-Bi satisfies ( 1441) . A non-degenerate 
holomorphic matrix function having Monodromy Properties (144D must be proportional to I. 
Since variables s m , s m , s m , s m are defined only up to a scalar factor, we can always choose 
them to make L = P^Ai-i □ 



The next theorem establishes a correspondence between Poisson brackets for Additive 
Representation fj43l) and those for a multiplicative representation. A proof using Krichever- 
Phong's Formula ([3]) is possible. However, we will take a different approach, since it provides 
a shorter proof. 
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Theorem 6. Let us define Poisson brackets on the direct product of d copies of single-pole 
spaces using Formulas (3l\): 

l^m' s mi "OO^m^mi \°r 

f § l § 3 | = ~2 r~2 
rm) ^mi ^m^mi rm 



|S° S 1 \ - -0 4 s 2 S 
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where we identify: 
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(46) 

the variables A m = u m+ i — u m and u m are constants from Theorem [3 Other brackets 
{s l m , s 3 m ,} vanish for m ^ m! . 

Then these brackets coincide with Quadratic Brackets (0|): 

{L{u) f L(u)} = [r(u - v), L(u) <g> L(v)}, 

where r(z) is elliptic r-matrix and L(z) = B d B d _i...B 1 . Consequently, they also coincide 
with Krichever- Phong's Formula (TJ|) (n = 2) on symplectic leaves. 



Symplectic leaves on the direct product are determined by 2d conditions 8((s 
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+ (*m) )) = 0; t/ieir dimension is 2d. 



Proof. Although the formula for B m (z) does not coincide with (136]) after Identification fj46|) . 
one can check that 

f B m (v)} = [r(u-v), B m {u) ® B m (v)], 

where r(z) is elliptic r-matrix (|4]). It is a direct computation using six Riemann's theta 
relations: 
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id six others, obtained by exchanging 
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Consequently: 

{Bi(u) f £»} = <%[r(u - v), Bi(u) ® Bj(v)\, 
and the group property of the quadratic bracket ([16J) implies: 

{L(u) f L(u)} = [r(u - f ), L{u) <g> L(t»)]. 

Proposition 3.3 in [TT] and Formula (|4ip imply that these brackets coincide with a two-form 
given by Krichever- Phong's universal formula on symplectic leaves. □ 
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